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Abstract. The sphere packing problem asks for the greatest density of a 
packing of congruent balls in Euclidean space. The current best upper bound 
in all sufHciently high dimensions is due to Kabatiansky and Levenshtein in 
1978. We revisit their argument and improve their bound by a constant factor 
using a simple geometric argument, and we extend the argument to packings 
in hyperbolic space, for which it gives an exponential improvement over the 
previously known bounds. Additionally, we show that the Cohn-Elkies linear 
programming bound is always at least as strong as the Kabatiansky-Levenshtein 
bound; this result is analogous to Rodemich's theorem in coding theory. Finally, 
we develop hyperbolic lineaj: programming bounds and prove the analogue of 
Rodemich's theorem there as well. 



1. Introduction 

What is the densest arrangement of non-overlapping, eongruent balls in M" ? This 
problem has a long history and has been extensively studied [CS99], and it has 
strong connections with physics and information theory [CIO]. With the proof of 
Kepler's conjecture by Hales [H05], the sphere packing problem has been solved in 
up to three dimensions, but no proof of optimality is known in any higher dimension, 
and there are only a few dozen cases in which there are even plausible conjectures 
for the densest packing. In and R'^* there are upper bounds that are remarkably 
close to the densities of the Es and Leech lattices, respectively; for example, Cohn 
and Kumar [CK04, CK09] came within a factor of 1 + 10^^^ of the density of Es 
and a factor of 1 + 1.65 • lO^''*' of the density of the Leech lattice. However, in most 
dimensions we must be content with much cruder bounds. In this paper, we will 
slightly improve the best bounds known in high dimensions, show how to obtain 
them via linear programming bounds, and extend them to hyperbolic space. 

The density of a sphere packing in M" is the fraction of space covered by the 
balls in the packing. More precisely, let B^{x) denote the ball of radius R centered 
at x; then the density of a packing is the limit as i? — > cxd of the fraction of B^{x) 
covered by the packing (the limit is independent of x if it exists). Of course this 
limit need not exist, but one can replace it with the upper density defined with a 
limit superior, and it is not hard to show that the least upper bound of the upper 
densities of all sphere packings in M" is actually achieved as the density of a packing. 
Let Art. denote this maximal packing density. 

A spherical code in dimension n with minimum angle ^ is a set of points on the 
unit sphere in M" with the property that no two points subtend an angle less than 
9 at the origin. In other words, {x, y) < cos 9 for all pairs of distinct points x,y in 
the spherical code. Let A{n, 9) denote the greatest size of such a spherical code. 
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In this paper, wo consider the problem of finding upper bounds for packing 
density. Linear programming bounds have proven to be a powerful tool. This 
technique was first developed by Delsarte [D72] in the setting of error-correcting 
codes, and his method can be extended to many other settings. In particular, 
Delsarte, Goethals, and Seidel [DGS77] and Kabatiansky and Levenshtein [KL78] 
independently formulated a linear program for proving upper bounds on A{n,9). 
Using this approach, Kabatiansky and Levenshtein found excellent upper bounds 
on A{n, 6) for large n, and they then applied a geometric argument to deduce a 
bound on Ar™ . Their upper bound is currently the best bound known for n > 43 
(see [CS99, p. 20]). It has the asymptotic form 

(1.1) Ar. < 2-(0-5990-+°(l))". 

Cohn and Elkies [CEOS] found a more direct approach to bounding sphere packing 
densities, with no need to consider spherical codes. Their technique set new records 
in every case with n > 4 for which the calculations were carried out; in particular, 
it improved on the previous record by Rogers [R58] in the entire range 4 < n < 42, 
although only n < 36 was reported in [CE03]. (Sec also Theorem 1.4 in [L0V12] 
for further improvements in low dimensions.) However, despite the evidence from 
low dimensions, the asymptotic behavior of the Cohn-Elkies bound is far from 
obvious and it has been unclear whether it improves on, or even matches, the 
Kabatiansky-Levenshtein bound asymptotically. Until this paper, it was only known 
how to use the Cohn-Elkies linear program to match the "second-best bound" by 
Levenshtein [L79] (sec Section 6 of [CE03]). 

The purpose of this paper is fourfold. In Section 2 we improve the Kabatiansky- 
Levenshtein bound by a constant factor by giving a simple modification of their 
geometric argument relating spherical codes to sphere packings. (This does not 
change the exponential decay rate in bound (1.1)). In Section 3 we show that in every 
dimension n, the Cohn-Elkies linear program can always match the Kabatiansky- 
Levenshtein approach. This further demonstrates the power of the linear pro- 
gramming bound for sphere packing. In Section 4 we prove an analogue of the 
Kabatiansky-Levenshtein bound in hyperbolic space. The resulting bound behaves 
the same as (1.1) asymptotically, and it is exponentially better than the best bound 
previously known in hyperbolic space. Finally, in Section 5, we develop the theory of 
hyperbolic linear programming bounds (based partly on unpublished work of Cohn, 
Luric, and Sarnak) and prove that they too subsume the Kabatiansky-Levenshtein 
approach. 

2. Geometric argument 

In all sufficiently high dimensions, the best upper bound currently known for 
sphere packing density is given by Kabatiansky and Levenshtein [KL78] (see also 
Chapter 9 of [CS99] and Chapter 8 of [Z99]). They first obtain an upper bound on 
A{n, 9) using linear programming and then use the inequality 

(2.1) Ar« < sin"(6i/2)A(n + 1, 9). 

The inequality was derived using a simple geometric argument. Here we improve it 
using an equally simple argument. 

Proposition 2.1. For alln>l and 60° < 9 < 180°, 

(2.2) ARn < sin"(6i/2)7l(n,6i). 
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Figure 1. Proof of Proposition 2.1. 



Since the unit sphere in M" can be embedded in the unit sphere in K"+-'^ via a 
hyperplanc through the origin, we always have A{n,9) < A{n + 1,9), with strict 
incquahty when 9 < 90°. The applications of (2.1) have 60° < 9 < 90°, so 
Proposition 2.1 will be a strict (though small) improvement. Neither inequality is 
useful in low dimensions; for example, when n = 2 and 9 = 60°, Proposition 2.1 
shows that Ar2 < 3/2. However, these inequalities arc valuable in high dimensions. 

For the sake of comparison, let us first recall the proof of (2.1). 

Proof of (2.1). Suppose we have a sphere packing in K" of density A using unit 
spheres. Consider a sphere 5g in ]R"+^ of radius R (to be chosen later), and place 

the sphere packing in M" onto a hyperplanc through the center of S'^, with the 
packing translated so that at least Ai?" of the sphere centers are contained in 
S^. This is always possible by an averaging argument. Project the sphere centers 
onto the upper hemisphere of S^, orthogonally to the hyperplanc. The projections 
of the sphere centers are still at least distance two apart, and thus separated by 
angles of at least 9, where sm{9/2) = 1/R. Therefore, Ai?" < A{n +1,6*), which is 
the bound that we wanted to prove, and we can achieve any angle by choosing R 
accordingly. □ 

Our motivation for revisiting this argument is that it feels somewhat unnatural 
to lift to a higher dimension in the process. Our proposition shows that a stronger 
inequality can be obtained without going to a higher dimension. The proof is similar 
to the techniques of [HSTIO] and [BM07], but this application appears to be new. 

Proof of Proposition 2.1. See Figure 1. Suppose we have a packing of unit spheres 
in with density A. Let "S*^"^ be a sphere in M" of radius R < 2 (to be chosen 
later), located so that it contains at least Ai?" of the centers of the spheres in the 
packing but its center is not one of them. Such a location always exists, by an 
averaging argument: a randomly chosen location will contain an average of Ai?" 
sphere centers. Now, project the sphere centers from the packing onto the surface of 
S^~^ using rays starting from the center of S^^^^. It follows from the lemma below 
that the projections are separated by angles of at least 9, where sin(0/2) = 1/R. 
Therefore, Ai?" < A{n, 9), as desired, and we can achieve any angle of 60° or more 
using i? < 2. □ 

Note that the proof breaks down if i? > 2, because two projected sphere centers 

can even coincide. 



Lemma 2.2. Suppose i? < 2. If XYZ is a triangle with \XY\ > 2, \XZ\ < R, 
\YZ\ < R, then ZXZY > 9, where sin(6»/2) = 1/i?. 
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Figure 2. Proof of Lemma 2.2. 

Proof. First we give a pictorial proof, which is more illuminating. In Figure 2, the 
bounds \XZ\ < R and \YZ\ < R place Z in the dark gray region, which is the 
intersection of the two disks centered at X and Y with radius R < 2. The light gray 
region contains all points P with ZXPY > 9. Since the dark region is contained 
inside the light region, it follows that ZXZY > 9. 

Let us also give an algebraic proof. Let x = \XZ\, y = \YZ\, z — \XY\, and 
7 = ZXZY. By the law of cosines, cos 7 = {x"^ + — z'^)/{2xy). By taking 
partial derivatives, we see that the expression (x^ + — z^)/{2xy) is maximized 
in the domain < x,y < R and z > 2 at (.t, y, z) = (R, R, 2). Therefore, cos 7 < 
1 - 2R-^ = 1-2 iim^{e/2) = cos 6. It follows that -f>6. □ 

Inequalities (2.1) and (2.2) can be stated a little more naturally in terms of 
packing density on the sphere. A spherical code on S"~^ with minimal angle 9 
corresponds to a packing with spherical caps of radius 9/2 with density 

, . . . fn sin""^ a; da; 

(2.3) A{n,e) \^ ■ n-2 , ■ 

sm X ax 

In other words, it covers this fraction of the sphere. Now if we let Agn-i(9) denote 
the optimal packing density, then (2.1) implies 

(2.4) -logAr < -logAs„-i(^^), 

n n 

where f{n) < g{n) mean limsup„_j.o2 f{n)lg{n) < 1. (At this level of asymptotics, 
the difference between (2.1) and (2.2) vanishes.) This simply amounts to verifying 
that 

1 /n/'sin"-2xdx 9 

- log % . n-2 1— ~ logsm -. 

n sm" ^xdx 2 

Furthermore, it is known that 

(2.5) -logAs„-l(^^) < -logAs„-.(</)) 
n n 

ioi 9 < (p (see (17) in [L75]). Thus, the exponential rate of the packing density 
for spherical caps is weakly increasing as a function of angle, and Euclidean space 
natiirally occurs as the zero angle limit. 

The proof of the Kabatiansky-Levenshtein bound (1.1) on Aa^i uses the following 
bound on A{n, 9) for < 9 < 90°, which is derived using the linear programming 
bound on spherical codes: 

,^ 1, ^, ,l+sin0, 1 + sin^ 1 — sin0, l — sm9 

2.6 - log A{n, 9) < . log — — log . . 

n 2sm6' 2sm6' 2sm^ 2sm^ 
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The bound (1.1) is then deduced by setting (2.6) into (2.1) and choosing 6 to 
minimize the bound, ^ which turns out to be around 9 « 63°. If we now apply our 
new inequahty (2.2) in place of (2.1), then we obtain an improvement in the bound 
by a factor of A{n+1, 6)/A{n, 6). If we use (2.6), then our bound is better by about 
a factor of 1.51 for the optimal value of 6. 

3. Linear programming bounds 

In [KL78] the upper bound on the maximum sphere packing density Ar™ was 
derived by first giving an upper bound for the maximum size A{n, 9) of a spherical 
code using hnear programming, and then using (2.1) to compare the two quantities. 
We refer to this method as the Kabatiansky-Levenshtein approach. Cohn and Elkies 
[CEOS] took a more direct approach to bounding Ar^, by setting up a different 
linear program. In this section, we show that the Cohn-Elkies linear program can 
always prove at least as strong a bound on Arti as the Kabatiansky-Levenshtein 
approach. 

This theorem is the continuous analogue of a theorem of Rodemich [R80] in 
coding theory (see Theorem 3.5 of [D94] for a proof of Rodemich's theorem, since 
Rodemich published only an abstract). Let A{n,d) denote the maximum size of a 
binary error-correcting code of block length n and minimal Hamming distance d 
(i.e., a subset of {0, 1}" with every two elements differing in at least d positions), 
and let A{n, d, w) denote the maximum size of such a code with constant weight w 
(i.e., every element of the subset has exactly w ones). The current best bounds on 
A(n, d) and A{n, rf, w) for large n are by McEliece, Rodemich, Rumsey, and Welch 
[MRRW77], using linear programming bounds. As in the Kabatiansky-Levenshtein 
approach, some of the best bounds on A{n,d) were obtained using bounds on 
A{n, d, w) along with an analogue of Proposition 2.1 known as the Bassalygo-Elias 
inequality [B65]: 

2^ 

(3.1) A{n,d) < j:^A{n,d,iu). 

The proof of (3.1) is by an easy averaging argument. In analogy with sphere packing, 
error-correcting codes play the role of sphere packings while constant weight codes 
play the role of spherical codes. Rodemich proved that any upper bound on A{n, d) 
obtained using the linear programming bound for A{n,d,w) combined with (3.1) 
can be obtained directly via the linear programming bound for A{n, d). Theorem 3.4 
below is the continuous analogue of Rodemich's theorem. 

3.1. LP bounds for spherical codes. We begin by reviewing the linear program- 
ming bounds for spherical codes. We follow the approach of Kabatiansky and 
Levenshtein [KL78], based on the inequality of the mean. 

Let 5""^ denote the unit sphere in M". A function /: [—1, 1] — > M is positive 

definite if for all N and all xi, . . . ,xn G S"'~^ , the matrix (/((x^, Xj))) ^^^^ ^.^^^ is 
positive semidefinite. (Note that this property depends on the choice of n; when 

^Lct us clarify a potentially confusing point. The fact that 6 63° minimizes the bound may, 
at first, seem to be at odds with (2.4) and (2.5), where we said that the exponential rate of the 
packing density Ay,,-i(6) is weakly increasing in 6. Both statements are correct. The bound in 
(2.6) is a preliminary bound on A{n,0), which can be improved for 6 less than the critical value 
(roughly 63°) by incorporating (2.5). 
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necessary for clarity, we will say such a function is positive definite on S" 
Equivalently, for all Xi, . . . , xn G 5""^ and ti,. . . ,tN € IR, 

X] titjfi{xi,Xj)) >0. 

l<i,j<k 

A result of Schoenberg [S42] characterizes continuous positive-definite functions as 
the nonnegative linear combinations of the Gegenbauer polynomials C^~^^^'^ for 

fc = 0, 1, 2, Recall that the polynomials C" are orthogonal with respect to the 

measure (1 — t"^)" dt on [—1, 1]. When a = (n — 3)/2, this measure arises naturally 
(up to scaling) as the orthogonal projection of the surface measure from S^~^ onto 
a coordinate axis. 

Given a positive-definite function g, define g to be its average 

__J^^g{t)il-f)(-^-^y^dt 

/!l(l-i2)("-3)/2rff 

with respect to this measure. Equivalently, g is the expectation of g{{x,y)) with x 
and y chosen independently and uniformly at random from S"*"^. If 

9{t) = Y.<^,ci'^~"^/\t), 

k>0 

then g = cq. 

Theorem 3.1 (Delsarte-Goethals-Seidel [DGS77], Kabatiansky-Levenshtein [KL78]). 
If g: [—1,1] — >■ K JS continuous and positive definite on 5""^, g{t) < for all 
l,cos^], andg>0, then 

9 

Let A^^{n, 9) denote the best upper bound on A{n, 9) that could be derived using 
Theorem 3.1. In other words, it is the infimum of g{'i-)/g over all valid auxiliary 
functions g. 

We will give a proof of this theorem following the approach of [KL78] , as prepa- 
ration for giving a new proof of Theorem 3.3 below. 

Proof. Let C be any spherical code in S""^ with minimal angle at least 9, let n be 
the surface measure on S"'^^, normalized to have total measure 1, let 6^ be a delta 
function at the point x, and let 

xec 

where A is a constant to be determined. We have 

jj g{{x,y))diy{x)du{y)>0, 

because we can approximate the integral with a sum and use the positive definiteness 
of g. This inequality amounts to 

X^g + 2X\C\g+ ^ g{{x,y))>0. 

x,yeC 
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Because (x, y) < cos 6 for distinct points x,y & C and g{t) < for ^ e [—1, cos 9], we 
have 

E 9{{x,y))<Y,9{{x,x)) = \C\g{l). 
x,yec xec 

Thus, 

X''g + 2X\C\g+\C\g{l)>0. 
To derive the best bound on \C\, we take A = —\C\. Then 

0<-|C|25 + |C|5(l) 

and hence 

\C\ < '-f, 

as desired. □ 

3.2. LP bounds in Euclidean space. The Kabatiansky-Levenshtein approach 
gives the following bound on Agn. The original version uses (2.1), but here we state 
the improved version using Proposition 2.1. 

Corollary 3.2. Suppose g satisfies the hypotheses of Theorem 3.1 with 60° < 9 < 
180°. Then 

Am„ < sin" (0/2) 4^. 

9 

Let us recall the Cohn-Elkies linear programming bound. Given an integrable 
function / : M" — > M, let / denote its Fourier transform, normalized by 



m = [ /(a;)e2-<-'*) dt. 



Let B'^ denote the n-dimensional ball with radius R. The volume of the n- 
dimensional unit ball is vol(Bf ) = 7r"/2/(n/2)!, where (n/2)! = r(n/2 + 1) for 
n odd. 

Much like the case of spheres, a function / : M" — >■ M is positive definite if for all 

N and all xi, . . . , x^ £ K", the matrix (/(x,; — Xj)'j is positive scmidefinite. 

A result of Bochner [B41] characterizes continuous positive-definite functions as the 
Fourier transforms of finite Borel measures. If / and / are both integrable, then / 
is positive definite if and only if / is nonnegative everywhere, by Fourier inversion 
and Bochner's theorem. 

Theorem 3.3 (Cohn-Elkies [CE03]). Suppose f : — > K is continuous, positive 
definite, and integrable, f{x) < for all \x\ > 2, and /(O) > 0. Then 

AMn<vol(Br)j|j- 

The original version in [CEOS] required suitable decay of / and / at infinity, 
and it was based on Poisson summation. These more restrictive hypotheses were 
removed in Section 9 of [CK07]. Here we give a more direct proof, although it has 
the disadvantage of not telling as much about what happens when equality holds as 
the Poisson summation proof does. 
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Proof. Without loss of generality, we can symmetrize to assume / is an even function 
(indeed, radially symmetric). This is not necessary for the proof, but it will simplify 
some of the expressions below. 

Let V he a packing with balls of radius 1, such that P has density Ar^. Given a 
radius r > 0, let Sr be the set of sphere centers from V that lie within the ball of 
radius r about the origin, let Vr be the volume of that ball, and let Nr = \Sr\- Then 

N 

lim voKBHTf = Am". 

Let R — r + ^ (in fact, -^/r could be replaced with any function that tends to 
infinity but is o(r)). Consider the signed measure 

where is the delta function at x, jij^ is Lebesgue measure on the ball of radius 
R centered at the origin, and A is a constant to be determined. As in the proof of 
Theorem 3.1, 

because / is positive definite. Equivalently, 
A'// f{x-y)dxdy + 2Xy2 [ f{x-y)dy+ V f{x-y)>0. 

Because f{x — y)<0 whenever x and y are distinct points in the packing, 

A'// f{x-y)dxdy + 2Xy^ [ f{x-y)dy + Nrf{0)>0. 

JJ\^\,\v\<R ^es^^l!/|<fl 

We set A = —Nr/Vr and divide by Nr to obtain 

Nr 1 /•/• ... ^ N 

Vr ' Vr 



II f{x-y)dxdy-2.^~y2 I f{x-y)dy + m>0. 

JJ\x\,\y\<R Vr '^rf-^J\y\<R 



It is not hard to compute the limits 

lim ^11 f{x-y)dxdy = m 

r^oo Vr JJ\x\,\v\<R 



^ E / f{x-y)dy = m. 

J\v\<R 



and 

lim 

Specifically, when < r, the y-integral covers all values of a; — y up to radius 
R — r = \fr. As r — )• DO these y-integrals converge to /(O), and all but a negligible 

fraction of the values of x satisfying \x\ < R also satisfy \x\ < r. 
Thus, in the limit as r — > oo we find that 

which is equivalent to the desired inequality. □ 
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Let denote the optimal upper bound on Ar™ using Theorem 3.3. Recall that 
A^^{n, 9) denotes the optimal upper bound on A{n, 9) obtained using Theorem 3.1. 
Our next result compares the LP bound on the sphere packing density Art, obtained 
from Corollary 3.2 with the one from Theorem 3.3. 

Theorem 3.4. For 60° < 9 < 180° and positive integers n, 

A^P < sm"{9/2)A^^{n,9). 

To prove Theorem 3.4, we will show that for any upper bound on Agn obtained 
using a function g in Corollary 3.2, we can always find a function / that gives a 
matching bound using Theorem 3.3. In other words, 

sin"(^/2)^=vol(Br)£M. 

9 /(O) 

We have a similar conclusion for the original Kabatiansky-Levenshtein bound 
using (2.1) without the 6 > 60° assumption. See the remarks following the proof. 

Proof of Theorem 3.4. Let g be any function satisfying the hypotheses of Theo- 
rem 3.1. The idea is to construct a function /: M" — >■ M based on g mimicking the 
geometric argument in the proof of Proposition 2.1. Let R= 1/ sin(6'/2), as in that 
proof. 

Consider the integral 



lB^{x)nB^{y)\{x,v} \\\X-Z\ W - Z\ 

where B'^{x) is the ball of radius R centered at x. Note that 

X- z y- z \ 

I 1 , I 1 ) = cos Zxzy, 

\x-z\ \y-z\l 

where Z.xzy denotes the angle at z formed by x and y. This angle is not defined if 
X = z OY y = z., but these cases occur with measure zero. 

The integral depends only on |a; — y|, so there is a radial function /: R" ^ M 
satisfying 

J Bl(^x)r^B^^{y)\{xM \V-AI ) 

We claim that / is positive definite. Indeed, let Xfl denote the characteristic function 
of S^(0) \ {0}. Then we can rewrite / as 

For any xi, . . . , Xjv G K" and ti, . . . ,tN G M, we can expand 

l<i,j<N 

as 

This expression is nonnegative, because g is positive definite on the unit sphere in 
R" and we can use tiXnixt — z) as coefficients. This shows that / is positive definite 
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on M". It is also iiitegrable, because it has compact support (it vanishes past radius 
2R). 

If \x\ > 2, then by Lemma 2.2, 

^ " < cos 6' 



\x-z\ \y-z\ 

for all z e n ^^(y) \ {x, y}. Since g{t) < for all t € [-1, cos 0] by hypothesis, 

it follows that f(x) < whenever \x\ > 2. Thus, / satisfies the hypotheses of 
Theorem 3.3. Now, /(O) = Yo\{B^)g{l) and 

/(0)= / f{x-0)dx 
= yol{B^fg. 

Therefore, 

vol(i.r)5^ = voKBD^^^ = =sin"(./2)^, 

as desired. □ 

When 9 < 60° wc can similarly match the Kabatiansky-Levenshtein bound 
obtained using (2.1) by adapting the above proof for the corresponding geometric 
argument. Let tt: _B" M"+^ denote the map that orthogonally projects the unit 
disk in the hypcrplanc in M"+^ to the upper half of the unit sphere in For 
any g in Theorem 3.1 that gives a bound for A{n + 1, 0), let 

/(x-»)=/ '(('(^l.'C^)))"- 

A similar argument shows that / is positive definite and f{x) < whenever |a;| > 2. 

We have /(O) = vol(Sg)g(l) and /(O) = vol(B]^)2E[g((7r(w), 7r(t>)))] > Yol{B'^)^g, 
where u and v are independent uniform random points in and the last inequality 
follows from g being positive definite on the unit sphere in ]R"+^ (using the inequality 
of the mean). However, wc cannot conclude that f{0) = vol{B'^)^g, so the version 
of this argument in Theorem 3.4 is more elegant. 

4. Hyperbolic sphere packing 

Hyperbolic sphere packing is far more subtle than Euclidean sphere packing. 
In both hyperbolic and Euclidean spaces, one must deal with the infinite volume 
of space available. The Euclidean solution is fairly straightforward: restrict to 
a large but bounded region, and then let the size of this region tend to infinity. 
The boundary eSiects have negligible influence on the global density. However, 
these arguments become much trickier in hyperbolic space, since the exponential 
volume growth means the limiting behavior is dominated by what happens near 
the boundary. Troubling phenomena occur, such as packings that have different 
densities when one uses regions centered at difii'erent points. There are numerous 
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other pathological examples (see, for example, Section 1 of [BR04]), and it is only 
recently that a widely accepted definition of density has been proposed by Bowen 
and Radin [BROS, BR04]. Before this definition, some density bounds were proved 
using local arguments that would apply to any reasonable definition of density, and 
indeed they apply to the Bowen-Radin definition. 

The best bound known is due to Boroczky [B78], who gave an upper bound for 
the fraction of each Voronoi cell that could be covered in a hyperbolic sphere packing. 
The bound depends on the radius of the spheres in the packing (the curvature of 
hyperbolic space sets a distance scale, so density is no longer scaling-invariant, 
as it is in Euclidean space). At least in sufficiently high dimensions, the bound 
is an increasing function of radius [M99] , so this bound is never better than the 
radius-zero limit. In that limit it degenerates to the Rogers bound [R58], which is 
asymptotic to 2~"/^ • n/e in n dimensions. 

Here, we improve the density bound to the Kabatiansky-Lcveiishtcin bound, 
regardless of the radius. Let Ah" (r) denote the optimal packing density for balls 
of radius r in H". We can bound the packing density of balls in hyperbolic space 
by the packing density of spherical caps on a sphere, as in the Euclidean setting 
discussed in Section 2. The next result is analogous to Proposition 2.1. 

Theorem 4.1. For all n > 1, 60° < < 180°, and r > 0, we have 

AH"(r) <sin"-i(0/2)A(n,(?). 

More precisely, one could replace sin"~^(^/2) with the hyperbolic volume ratio 
vo1(B;?)/vo1(S^), where R is defined by sinhi? = (sinhr)/ sin(e/2). That would 
slightly improve the inequality without changing the proof, at the cost of making 
the statement more cumbersome. 

As in the Euclidean case (2.4), this theorem implies that 

(4.1) sup-logAH"(r) < -logAs„-i(^). 

r>o n n 

By using the Kabatiansky-Levenshtein bound on Agn-i, i.e., (2.6) with « 63°, we 
obtain the following new bound on Ajjn , which is an exponential improvement over 
the previous best bound due to Boroczky [B78]. It is independent of the radius of 
the balls used in the packing. 

Corollary 4.2. We have 

supAHn(r)<2-(o-5990-+<'(i))». 

r>0 

The Bowen-Radin approach to hyperbolic packing is based on ergodic theory, 
but our argument is elementary. All we need is the following fact: if there exists a 
packing of density 6 using balls of radius r, then for every i? > 0, there exists a ball 
B of radius R containing a subset of at least 

^vol(B|) 
vol(B») 

points at distance 2r or more from each other and not equal to the center of B. 
Naively, this should follow from a simple averaging argument, since if we place B at 
random in a dense packing, then this is the expected number of sphere centers it 
will contain, and the probability that one of them will hit the center of B is zero. In 
Appendix A, we explain why the Bowen-Radin density has this property. 
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The proof of Theorem 4.1 is analogous to the Euchdean case. The heart of the 
proof is the following lemma. 

Lemma 4.3. Let r < R < 2r and sin | = ^^'^ • a packing of spheres of radius 
r m H", every ball of radius R contains at most A{n,9) sphere centers other than 
its own center. 

Proof. We use the same projection argument as in the proof of Proposition 2.1. 
Project the sphere centers from the packing onto the surface of the ball of radius R 
using rays starting from the center of the ball. By the next lemma, the projections 
are separated by angles of at least 6, so there can be at most A{n, 9) of them. □ 

The next lemma is the hyperbolic analogue of Lemma 2.2. 

Lemma 4.4. Consider a hyperbolic triangle with side lengths a, b, c, and the angle 
opposite to c having measure 7. If < a,b < R < 2r < c, then 

7 sinh r 
sin - > — — . 
2 ~ sinhi? 

Proof. By hyperbolic law of cosines, 

cosh a cosh b — cosh c 
sinh a sinh b 

Let 

, , cosh a cosh 5 — cosh c 

f[a,b,c) = — r-TT • 

smh a smh 

We wish to maximize f{a, b, c) in the domain < a,b < R < 2r < c. Since / is 
monotonically decreasing in c, it is maximized by setting c = 2r. We have 

df cosh a cosh c — cosh b 

da sinh^ a sinh b 

which is nonnegative since coshc > cosh 6 and cosh a > 1. Thus f{a,b,c) is 
nondccrcasing in a, and it is maximized by setting a = R. The same is true for b by 
symmetry, and so 

, . cosh^ i? — cosh 2r ^ 2sinh'^r 

cos 7 = / (a, 0, c) < 2 = 1 2 — • 

sinh R sinh R 

Therefore 

. n 7 1 — COS 7 sinh^ r 

2 2 "sinh^ii' 

and the result follows. □ 

Proof of Theorem 4.1. Define R to satisfy sin | = Since 60° < 61 < 180°, we 

have r < R<2r. (Note the latter inequality does not always hold when 9 < 60°. It 
fails in the limit as r ^ but holds for large r.) 

Given a packing of balls of radius r in M" with density S, there exists a ball B of 
radius R containing at least sphere centers from the packing that are not 

the center of B itself. By Lemma 4.3, 

-^^<-'^ 

The volume of a ball in H" is given by 

vo1(B;?) = f2„ / sinh"-^ x dx, 
Jo 
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where fin = 27r"/^/r(n/2) is the surface volume of the unit Euchdean (n — l)-sphere. 
Thus 

, , rr sinh"'""'^ xdx , 

(4.2) 

/ sinhr 



ra-l 



- Vsinhi?y ^ ' ' 

= sin"-i(6l/2)^(n,6l), 
wlicrc the second incquahty follows from Lemma 4.5 below. □ 

If we fix the ratio ^3^^]^ , then the ratio of the integrals in (4.2) is almost determined 
by the following lemma (the lower bound is sharp as r — )• and the upper bound is 
sharp as r — > oo). Wc do not need the lower bound, but it shows that (4.1) cannot 
be substantially improved by a more careful analysis of the volume of hyperbolic 
balls. 



Lemma 4.5. For < r < R, 



( sinh r \ " \l sinh" ^ xdx ( sinh r ^ " ^ 



V sinh R ) sinh"- ^xdx \ sinh R 

Proof. These inequalities amount to saying that 

sinh"-^ a;rfa; 
sinh"-^ r 

is an increasing function of r, while 

/J^sinh"^^ xdx 
sinh" r 

is a decreasing function of r. 

The derivative of the former function is 

(n - 1) coshr , , ^ 

1 —j: / smh X dx, 

smh r 

so we must prove that 

' sinh"-ixdx>0. 



/ 

Jo 



(n — 1) coshr 

The left side of this inequality vanishes when r = 0, and its derivative with respect 
to r is 

sinh"-V 
(n — 1) cosh^ r ' 
so it is increasing and hence nonnegative. 
To show that 

Jq sinh" a; da; 
sinh" r 

is decreasing, note that its derivative is 

1 n cosh r 



sinh r sinh"+^ ' 



/'sinh"-i 
Jo 



xdx, 
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SO we must prove that 

n cosh r 

Again the left side vanishes when r = 0, and this time its derivative is 

sinh"+^ r 
n cosh^ r ' 

so it is increasing and hence nonnegative. This completes the proof. □ 

5. Linear programming bounds in hyperbolic space 

It is natural to try to extend the results of Section 3 on linear programming 

bounds to hyperbolic space, but one runs into technical difficulties. 

Given a function / : [0, oo) — K, we view it as a radial function on hyperbolic 
space and define the corresponding kernel /: H" x H" ^ M by f{x,y) = f{d{x,y)), 
where d denotes the metric on H". (Using the same symbol for both functions 
is an abuse of notation, but it is convenient not to have to write the metric d 
repeatedly, and the number of arguments makes it unambiguous.) We say / is 
positive definite if for all N and all Xi, . . . ,xpf E H", the matrix [f{xi, Xj)) 

is positive semidefinite, and we say it is integrable on H" if a; f{x,y) is an 
integrable function on H" (of course this is independent of y). 

Let G be the connected component of the identity in the isometry group of H", 
and let K be the stabilizer within G of a point e e H". Then (G, K) is a Gelfand 
pair; i.e., the algebra L^{K\G/K) of integrable, bi-Zi'-invariant functions on G 
forms a commutative algebra under convolution. Here G/K is H" and functions 
on K\G/K correspond to radial functions on H". See Chapters 8 and 9 of [W07] 
for an account of Gelfand pairs and spherical transforms (and see [T82] for a more 
concrete exposition of Fourier analysis in H^). In the setting of H", this theory 
gives a well-behaved Fourier transform for radial functions. For each A > 0, let Px 
be the unique radial eigenfunction of the Laplacian on H" with eigenvalue A and 
^'a(O) = 1. These functions arc positive definite for all A > (sec Theorem 5.2 in 
[T63, p. 346]). Given a function /: [0,oo) — >• M that is integrable on H", its radial 
Fourier transform is given by 

/(A)= / f{x,e)Px{x,e)dx, 

which is of course independent of e G H". As in the Euclidean case, the Fourier 
transform extends to L?(K\G/K), and it yields an isomorphism from that space to 
X^([0, oo),/ip), where hp is the Plancherel measure. However, unlike the Euclidean 
case, the Plancherel measure for H" is supported just on [n^/4,oo). 

Positive-definite functions are characterized by the Bochner-Godement theorem 
(see Theorems 9.3.4 and 9.4.1 in [W07] or Theorem 12.10 in Chapter III of [H08]). 
For continuous, integrable functions, it says that / is positive definite if and only 
if / is nonnegative on the support of the Plancherel measure. However, / can 
be negative outside of the support, because G is not amenable: Valette [V98] 
has constructed a continuous, positive-definite function with compact support and 
negative integral. (His construction works in G, rather than G/K, but it is easy to 
make it bi-ii'- invariant.) 
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In the linear programming bounds, we will assume / > everywhere, which is a 
strictly stronger assumption than positive definiteness. We do not know whether 
the stronger hypothesis is truly needed for the following conjecture, but it will be 
needed for the proof of Theorem 5.5. 

Conjecture 5.1. Let f: [0,00) — J- M 6e continuous and integrable on H", and 
suppose f{x) < for all x >2r while /(A) > for all X> and /(O) > 0. Then 

Au^{r)<vom&. 

/(O) 

Here, of course, vol(i?") denotes the volume of a ball of radius r in H". 

Let (r) be the best bound on Ah^ (r) obtainable by optimizing the choice of 
/ in Conjecture 5.1. The following theorem is the hyperbolic analogue of Rodemich's 
theorem. 

Theorem 5.2. For 60° < 9 < 180°, positive integers n, and r > 0, 

ALP(r)<sin"-i(0/2)ALP(n,0). 

Proof. The argument is much like the proof of Theorem 3.4. Define R by sinh R = 
(sinhr)/ sin(0/2), and let XR denote the characteristic function of a ball of radius 
R. Given a function g satisfying the hypotheses of Theorem 3.1, define / by 

f{x,y)= / XR{x,z)xR{y,z)g{cosZxzy)dz, 



where Zxzy denotes the angle at z formed by the geodesies to x and y. Of course, 
this angle is not defined when x = z oi y = z, but these cases occur with measure 
zero. 

Exactly the same approach as in the proof of Theorem 3.4 shows that / is a 
positive-definite function and that f{x, y) <0 when d{x, y) > 2r. However, merely 
being positive definite does not imply that f{X) > for all A > 0. To prove that, 
we start by fixing y G H" and writing 

/(A)= / P^{x,y)f{x,y)dx 

Xr{x, z)xR{y, z)Px{x, y)g{cos Zxzy) dx dz. 



The integrand now depends only on d{x,z), d(i/,z), and Zxzy (note that the 
hyperbolic law of cosines determines d{x,y) using this data). If wc fix z instead of 
y and integrate over x and y, then the induced measure on these three parameters 
is the same. Thus, for fixed z we have 

fW= / XRix,z)xR{y,z)Px{x,y)g{cosZxzy)dxdy. 

Now wc can see that /(A) > 0, because {x,y) Px{x,y)g {cos Zxzy) defines a 
positive-definite kernel for x,y gW^\ {z}. Specifically, the product of two positive- 
definite kernels is positive definite by the Schur product theorem (Theorem 7.5.3 in 
[HJ13]), which says that the set of positive-semidefinite matrices is closed under the 
Hadamard product. 

It also follows from this formula that 

/(O) = Yol{B^)\ 



16 



HENRY COHN AND YUFEI ZHAO 



and combining this with /(O) = vol(_B^)5(l) yields 

/(O) vol(B^) g ^ ' ' g ' 

as desired. □ 

In the remainder of this section, we explain why a straightforward approach fails 
to prove Conjecture 5.1 and how to prove it for periodic packings. The latter proof 
is based on unpublished work of Cohn, Lurie, and Sarnak. 

5.1. Obstacles to proving the conjecture. Wc have not been able to prove 
Conjecture 5.1 by imitating the proof of Theorem 3.3. The problem is that the 
boundary effects when restricting to a ball are not negligible. 

Specifically, consider a hyperbolic sphere packing with balls of radius r, and 
imagine restricting it to a ball of radius R (i.e., looking only at the points within 
this large ball). Let Sr be the set of all sphere centers in this ball and the 
Lebesgue measure on the ball, and consider the signed measure 

where A is a constant to be specified shortly. Because / is positive definite, 



f{x,y)dv{x)dv{y) > 0, 

J J 

which implies 

JJ fix, y) dimix) duRiy) + ^>^T. J ^'"^(y) + 1^^1/(0) ^ 

xGSr 

as in the proof of Theorem 3.3. Under the Bowcn-Radin definition of density, saying 
there is a packing with density 6 means there is an isometry-invariant measure on 
such packings (see Appendix A for an overview of the Bowen-Radin theory). If we 
average over such a measure, we find that 

2SX 



'^)JJ fix,y) df,n{x)df,R{y) + S^^l^f {0)>0. 



vol(S. 

In particular, taking A = — (5/vol(i3") yields 

/(O) 



5 < voim ^^^^ ^^^^^^ ^^^^^^^ . 

This proves a legitimate bound on the density: 

Proposition 5.3. Let f : [0,oo) — > M be conUnuous and positive definite on H", 
and suppose f{x) < for all x > 2r. Then for each R> 0, 

AMr) < -l(^")jJ;(,,,)rf^^(fSX^(,)/,ol(Bn)' 

assuming the denominator is not zero. 

If is natural to add the assumption that / is integrable and then take the limit 
as — >■ oo. However, the denominator does not converge to /(O), as it does in the 
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Euclidean case. To see why, let XR be the characteristic function of a ball of radius 
R. Then 

/ / f{x, y) dfj.R{x) diiR{y) = (/ * Xr)Xr, 
J J Jm"- 

where the right side denotes the integral of a radial function on H". Because the 

Fourier transform is unitary, this integral equals 



f^XRXRdlJ-p, 

where hp is the Planchercl measure, and that simplifies to 

fx%dl^p. 

One can show similarly that vol(S^) = J Xr dup, and thus 

II f{x,y)dfj,R{x)dnR{y) ^ / fx\ djJ-p 
vol(BS) Ixldf^p ■ 

We can already sec a problem: we would like the mass of Xr to be concentrated near 
as i? — )■ oo. However, is not even contained within the support of the Plancherel 
measure /xp, so this cannot possibly work. To see how badly it fails, we return to 
radial functions on H" via 



/ fXR dnp= f ■ ixR * Xr)- 

J JM" 



(We use • for multiplication here to avoid confusion with / applied to an argument.) 
The function {xr * Xr) I vol(B^) measures the fraction of overlap between two balls 
of radius i? whose centers are a given distance apart. In the limit as i? — > oo, this 
function does not converge to 1, as it does in Euclidean space. Instead, when the 
distance between the centers is z it converges to 

D ( 1 . n— 1 n— 1 \ 

(5.1) ^ ' ' ^ 



r> / 1 . n— 1 n— 1 ^ 
V2' 2 ' 1 , 



where 

B{u\ a, 



Jq 



is the incomplete beta function. (See Appendix B for the calculation.) Note that 
(5.1) equals 1 when z = and vanishes in the limit as ^ — )• oo. 
Thus, 

II fix,y)dfiR{x)diJ.R{y) 

converges to the integral of / times a function that takes values between and 1. 

Variants of this approach, for example replacing xr with a smoother function such 
as the heat kernel, fail for essentially the same reason. In Euclidean space the heat 
kernel converges to the constant function 1 as time tends to infinity, if we rescale it 
so its value at the origin is fixed as 1. In other words, flowing heat becomes nearly 
uniformly distributed over time. However, in hyperbolic space that is not true (heat 
kernel asymptotics can be found in [DM88]). 
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5.2. Periodic packings. We can prove a variant of Conjecture 5.1 in the special 
case of periodic packings, i.e., packings that are invariant under a diserete, finite- 
covolume group of isometries. This was first proved by Cohn, Lurie, and Sarnak in 
unpublished work; here, we give a proof under weaker hypotheses but using the same 
fundamental approach. It is not known whether periodic packings come arbitrarily 
close to the Bowen-Radin packing density, although it is known that they cannot 
always achieve it exactly. Maximizing the density for a single-orbit packing in H" is 
equivalent to maximizing the systolic ratio of an n-dimensional hyperbolic manifold 
with finite volume (see [K07] for background on systolic geometry), which makes 
periodic packings a particularly important case. 

Given a periodic packing with balls of radius r, let F C G be its symmetry group. 
By assumption, r\H" has finite volume. Suppose the spheres in the packing are 
centered on the orbits ra;i, . . . , Txn, and let Fj be the stabilizer of Xi in F. Then the 
density of the packing is the fraction of a fundamental domain covered by balls. If 
the stabilizers are trivial, then this fraction is simply A^vol(i?")/ vol(F\H"). More 
generally, F^ preserves the ball centered at Xi, and only a l/|Fi| fraction of this ball 
will lie in any given fundamental domain (specifically, one element of each Fj-orbit). 
Thus, the density is 

voi(g;^) ^ 1 

vol(F\H") ^ |Fi|' 

It is the same as the density under the Bowen-Radin definition. 

The proof of the linear programming bounds will be based on the Selberg trace 
formula [S56] (or, more precisely, a pre-trace formula). Selberg did not publish 
a complete proof of his formula. For a detailed proof, see [F67] for H^, [V73] for 
hyperbolic spaces under some mild hypotheses on the discrete group, [EGM98] for 
the three-dimensional case (using techniques that work in greater generality), and 
[CS80] for hyperbolic spaces in general. 

Call a radial function / admissible if it is continuous and integrable on H" and 
for every discrete subgroup F C G for which F\H"' has finite volume, 

converges absolutely for all x,y G H" and uniformly on compact subsets of H". For 
example, by Lemma 1.4 in [V73] (attributed to Selberg), / is admissible if it is 
integrable on H" and there exist constants ci,C2 such that 

\f{x,y)\<ci [ \f{x,z)\dz. 

J d{z,y)<C2 

The only place where we require the trace formula machinery is the proof of the 
following lemma: 

Lemma 5.4. Let f: [0,oo) — > K 6e admissible on H" and satisfy /(A) > for all 
A > 0. If T is a discrete subgroup in G with finite covolume, then the function F 
defined by 

y) = ^ filx, y) 

is positive definite on H" x H". Furthermore, F — /(O)/ vol(F\]HI") remains positive 
definite. 
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In other words, for all xi, . . . , xat € H", the NxN matrix with entries F{xi, xj) — 
/(O)/ vol(r\]HI") is positive semidefinite. 

Proof. First, suppose r\]HI" is compact, so the spectrum of the Laplacian is discrete. 

Let vo,vi, . . . be the orthonormal eigenfunctions of the Laplacian on r\H", viewed 
as periodic functions on H", and let Aq < Ai < . . . be the corresponding eigenvalues. 
The sum 

is periodic modulo F as a function of x (or y), so we can expand it in terms of the 
eigenfunctions of the Laplacian. We have 



V/(7a;,y) ~ Vt>i(a;) / V. 
7er i=o -^rXHn ^^^^ 



',f{iz,y) I v,{z)dz, 
7er i=o \™" \7er 

where ~ denotes L"^ convergence. The coefficients unfold to 



f{z,y)vi{z)dz, 



and we can rotationally symmetrize about y, which turns Vi{z) into Vi{y)P\^{z,y) 
and the coefficient into 

Vi{y) / f{z,y)Pxi{z,y)dz = Vi{y)f{Xi). 
(The conjugate on P\. does not matter, because this function is real-valued.) Thus, 

oo 

^f{ix,y) c^^f{Xi)vi{x)vi{y). 

jeT i=0 



The functions (a;, y) Vi{x)vi{y) are clearly positive definite, and the coefiicients 
/(Ai) arc; nonncgative. Furthermore, positive dcfinitoncss is preserved under point- 
wise convergence. However, this expansion may not converge pointwise. Fortunately, 
L"^ convergence implies that a subsequence converges pointwise almost everywhere 
(see, for example, Theorem 3.12 in [R87]). Thus, for almost all .xi, . . . ,xm G H", 
the matrix with entries F{xi,Xj) is positive semidefinite, and the same holds for 
all a;i, . . . ,xn by continuity. Furthermore, vo is the constant eigenfunction, so vqVo 
must be l/vol(F\H") by orthonormality, and thus F - /(O)/ vol(r\H") is also 
positive definite. 

All that remains is to deal with the case when T\W" has finite volume but is not 
compact. Harmonic analysis on the quotient is quite a bit more involved, because of 
continuous spectrum coming from the cusps, but a completely analogous argument 
works. Suppose r\IHI" has h cusps. For 1 < k < h and s G C with 5R(s) = n/2, 
there is an Eisenstein series x i— >■ Ek{x, s), which is an eigenfunction with eigenvalue 
s{n — s). When s = n/2 + it, the eigenvalue becomes + 1"^, so it is contained in 
the support oo) of the Plancherel measure. The spectral resolution is now 



X /(7a;, y)c^Yl /(^j)^j (y) 
7er j=o 

" )dt 

fe=i" 



+ Y. /(nV4 + f)Ej {x, n/2 + it)Ej {y, n/2 + it) , 
1 1 J— 00 
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(see (7.30) in [CS80, p. 75]), and the proof precedes as in the compact case. □ 

The proof of Lemma 5.4 depends on the hypothesis that / > everywhere, not 
just on the support [n-^/4,oo) of the Plancherel measure, because the eigenvalues of 
r\]HI" need not be contained in the support: is always an eigenvalue, and there 
can be others between and Selberg's eigenvalue conjecture [S95] says there 

are no nonzero eigenvalues in (0, 1/4) for congruence subgroups of SL2(Z), but that 
is a special arithmetic property not shared by more general groups. 

Theorem 5.5 (Cohn, Lurie, and Sarnak). Let /: [0, oo) — >• M 6e admissible on H", 
and suppose f{x) < for all x >2r while /(A) > for all X > and /(O) > 0. 
Then every periodic packing in H" using balls of radius r has density at most 

vol(5?)m. 
/(O) 

Proof. Consider a periodic packing consisting of orbits r.xi, . . . .Tx^ of a finitc- 
covolume group F, and let Fj be the stabilizer of Xi in F. Recall that the density of 
the packing is 

N 

\ol(Bl'.) 



Let 



By Lemma 5.4, 



vol(F\H") ^ 

F{x,y) = ^f{ix,y). 
7er 



^3^jF,||r,| 1^^^^"^^^ vol(F\H"); - 



which amounts to 



^ 1 

X] rrTi 



/(o) 



|ril|r,-| - voi(r\H")- 

On the other hand, F{xi, xj) < for i ^ j, because all the terms in the sum defining 
F are nonnegative in that case, and F{xi,Xi) < |rj|/(0), because there are \Ti\ 
group elements 7 for which ^Xi = Xi. Thus, 

^|ri|/(o) ^ 1 V /(o) 



v iri|/(o) ^ V 1 
U - voi(r\H")- 



We conclude that 



as desired. □ 
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Appendix A. Center density in hyperbolic space 

In this appendix, we briefly sketch how our work fits into the Bowen-Radin 
framework [BROS, BR04]. Let Sr be the space of relatively dense packings of H" 
with closed balls of radius r (i.e., packings in which any additional such ball would 
intersect one from the packing). Bowen and Radin give a natural metric to Sr, under 
which it is compact, and they study the action of the isometry group G of H" on Sr- 
They define random packings by G-invariant Borel probability measures jj. on Sr, 
and they define the density of /j, to be the probability that the origin is contained in 
one of the balls in the packing. The optimal packing density Ah»i (r) is the least 
upper bound for the density of any such measure. Bowen and Radin show that this 
density is achieved by some measure, and they show how to obtain well-behaved 
dense sphere packings by sampling from such a distribution. They show that this 
definition encompasses the reasonable examples that were known before, such as 
packings with finite-covolume symmetry groups, and they make a convincing case 
that it is the right definition. See also [R04] for intuition and background. 

Let iJ. be an invariant probability measure on Sr, and suppose fi has density S. 
We would like to show that the expected number of sphere centers in a given ball of 
radius R is 

volOB|) 
vo1(B;?) ' 

as needed in Section 4. (Because jj, is G-invariant, the probability that the center of 
the ball of radius B will be one of them is zero.) By the definition of density, this is 
true for R = r. To prove it in general, we will use the Benjamini-Schramm mass 
transport principle, which is Theorem 5.2 in [BSOO]. That theorem is stated only 
for H^, but they explain on page 497 that it holds in general symmetric spaces. 

The mass transport principle says that if is a Borel measure on H" x H" that 
is invariant under the diagonal action of G, and if u{A x H") < oo for some open 
set A, then 

u{A X H") = u{ir X A) 

for all measurable A. 

Given a packing P G Sr, define a measure vp on H" x H" by 

X 

where x runs over all sphere centers in P, Xb"(x) denotes the characteristic function 
of the ball of radius r about x (viewed as a measure on H"), and 6^ is a delta 
function at x. Let v be the average of vp over all P G Sr, according to the measure 
H- Because fj, is G-invariant, so is v. 
Now 

z^(H" X ^) = vol(i?") • (expected number of sphere centers in A), 

while 

p{A X H") = vo\{A) ■ (expected number of sphere centers in a ball of radius r) 
= vol(A) • 5. 

By the mass transport principle, i/(A x H") = i/(H" x A), and taking A = 
completes the proof. 
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Appendix B. Overlap of balls in hyperbolic space 

Proposition B.l. If xi and X2 are points in H" at distance r from each other, 
then 



vol (Bg(xi) n B^{X2)) _ ^ (tT? 



n— 1 n— 1 

- 2 ' 2 

lim 



Recall that 

B{u;a,/3) = / f^-^il - tf-^ dt 
Jo 

is the incomplete beta function. To prove Proposition B.l, we will compute the 
convolution {xr * Xr)I vo1(-B]^) on H", where XR the characteristic function of a 
ball of radius i?, and take the limit as i? — > oo. 

First, we observe that given two radial functions fi and /2 on H", their convolution 
can be computed as follows. Let xi and X2 be two points in H" at distance r from 
each other, and consider a third point z at distances ri and r2 from xi and x^, 
respectively. Then (/i */2)(r) is the integral of fi{ri)f2{r2) over all z G H". To write 
it down explicitly, we can use polar coordinates centered at xi. Let u = cos Zx2Xiz. 
Then 

o (n-l)/2 poo pi 

(./i*./2)M = Y{(n^l)/2) X j Jiir,)f2{r2){smh--'r^)il-u')^--^^/Uudri, 

where we view r2 as a function of ri and u via the hyperbolic law of cosines. 
Changing variables from u to r2 yields 

/TiiN 27r("~^)/^ f f J. / Nj./ , sinhrisinhr2 „, ^(n-3)/2 i , 

where the integral is over all ri and r2 such that r, ri , r2 form the side lengths of a 
triangle, and 

C(r, ri, r2) = 1 — cosh^ r — cosh^ ri — cosh^ r2 + 2 cosh r cosh n coshr2. 
Note that C{r,ri,r2) can be factored as 

( cosh(ri + r2) — cosh r) ( cosh(r + ri) — cosh r2) ( cosh(r + r2) — cosh n) 

cosh(r + ri + r2) — 1 

Now we apply (B.l) to fi = f2 = xr divide by vol(_B]^), which is asymptotic 
to 27r"/2e("-i)^/((n-l)2"-ir(n/2)) as E ^ oo. We find that {xR*XR){r)l vol(Bg) 
is asymptotic to 

where P) = B(\;a, (3) = V{a)T(fi) /T{a + jS) denotes the beta function and X 
is the set of (ri, r2) S [0, i?]^ for which r, ri, r2 form a triangle. 

We now change variables in (B.2) from ri and r2 to a; = ri—r2 and y = 2R—ri—r2- 
In the new variables, the domain X of integration becomes 

{{x, y) : < r and |a:| < y < 2R — r}. 

As i? -S^ oo, we find that sinhn ~ e^+(^-2')/V2, sinhr2 - e^-i^+y)!-^ j^^ and 

, no „coshr — coshx 
C(r-,ri,r2) ~ e^^'^ ^ . 
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The integral (B.2) converges to 

(^-1)^'""'^^' f /"e-("-i)^/2^coshr-coshx)("-3)/2dydx. 

B((n-l)/2,l/2)sinh"-Vy_.yH ^ ' 

We can evaluate the y integral explicitly, and the remaining integrand is an even 
function of x. The integral thus equals 

9(n-l)/2 rr 

^ / e-("-iW2(coshr-cosha;)("-=^)/2da;. 

B((n- l)/2,l/2)sinh"-2r Jo 

Finally, we change to a new variable 

t = 



e' — e" 

to arrive at 



/ (i(l-t))^ 

Jo 



i3((n-l)/2,l/2) 7o 

It follows from the duplication formula for the gamma function that this expression 
equals 

B((n-l)/2,(n-l)/2)/2 ' 
This completes the proof of Proposition B.l. 
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